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For the implementation of a quantum computer it is nec- 
essary to exercise complete control over the Hamiltonian of 
the used physical system. For NMR quantum computing the 
effectively acting Hamiltonian can be manipulated via pulse 
sequences. Here we examine a register consisting of A'^ selec- 
tively addressable spins with pairwise coupling between each 
spin pair. We show that complete decoupling of the spins is 
possible independent of the particular form of the spin-spin 
interaction. The proposed method based on orthogonal ar- 
rays is efficient in the sense that the effort regarding time 
and amount of pulses increases only polynomially with the 
size A'^ of the register. However, the effect of external control 
errors in terms of inaccurate control pulses eventually limits 
the achievable precision. 
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I. INTRODUCTION 

Couplings within a quantum register may be consid- 
ered a necessary evil: the two-qubit operations needed 
for the realization of quantum computing ffl rely on in- 
teractions between different qubits. For fast gate opera- 
tions even strong couplings are desirable. On the other 
hand, residual interactions, that are always present, dis- 
turb the system and lead to unwanted evolution, so that 
stored information may get lost. This kind of "hard- 
ware" problems exist even in situations where external 
(classical) contacts are supposed to control pertinent pa- 
rameters of the network Hamiltonian |g,^ . They are par- 
ticularly relevant for conventional NMR implementations 
y,^ , where couplings cannot simply be turned on and off, 
but are predetermined by the respective Hamiltonian. 

The interplay between local disorder and interactions 
has recently been investigated with respect to the possi- 
ble emergence of the so-called "quantum chaos" [p| . As 
expected, disorder helps to avoid the delocalization of lo- 
cal register states. On the other hand, (static) disorder 
itself presents a kind of control loss, which is not wel- 
come. We will include some aspects thereof via faulty 
gate operations. 

Fortunately, by using sophisticated pulse sequences a 
wide range of control can be exerted on the system. Cou- 
pling terms can be removed in an NMR spin echo exper- 
iment by applying so-called refocusing pulses, essentially 
reversing the time evolution. This method can be ex- 
tended to the task of removing selected or even all cou- 
pling terms in a quantum register of N qubits. Such a 



scheme has been described by Linden et al. m. Unfortu- 
nately it lacks efficiency. Efficient schemes have been pro- 
posed by Leung et al. Q and Jones and Knill ||] , respec- 
tively. However, these studies focus on the special model 
of weak scalar coupling only. For other coupling models 
it remains open whether there exist generalized schemes, 
let alone efficient ones. In this paper we will demon- 
strate that, indeed, any interaction can be suppressed 
efficiently, regardless of the specific coupling model, and 
even without knowledge of the underlying Hamiltonian. 
Perfect decoupling in general requires fast application 
of many pulses. In a real implementation the pulses will 
never be ideal, which will impair the decoupling. Quan- 
tum error correction methods could be used, but they 
require supplementary resources in the form of auxiliary 
qubits. Here, we want to probe the limitations to the 
given method that arise from faulty gate operations with- 
out considering error correction methods. 



II. MODEL AND AVERAGE HAMILTONIAN 

We consider a system of N coupled spin- 1/2 particles 
(qubits) with the general Hamiltonian model 
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where the ai {i = 1,2,3 or x, y, z) denote the Pauli op- 
erators. The terms linear in those are the Zeeman terms 
specifiedby the 3-dimensional vectors a; ('^' = {uj] ). The 
interaction between two qubits /x and u is described by 
the couphng tensor H'^" — {H^"). We will call couplings 
"diagonal" when they are represented by a diagonal 
coupling tensor H = diag( Jj,, Jy, Jz). Examples are the 
strong and weak scalar coupling, with J^ — Jy = Jz —'■ J 
and Jx = Jy — 0, respectively. 

We will now subject the system to a cyclic pulse train. 
The pulses are assumed to be infinitely short, so that they 
can be represented by instantaneous operations i\. . Each 
pulse is followed by a free evolution period of duration t^ ■ 
A sequence consisting of n pulses will thus be denoted by 



{Pl,Ti,...,Pn,Tn), 



(2.2) 



where the pulses are applied from left to right. For 
the further analysis we use the formalism of average 
Hamiltonian (AH) theory |10|-[l^. The time evolution 



operator is given by U{tc) 
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"'^^"Pi., where 



Efc=i '''k i^ the cycle time (all products are ordered 



with increasing k from right to left). For a cychc se- 
quence (defined by nfe=i -Pfc — 1)' ^^ <^^^ write this as 

Uitc) = U'UmliPi)-'^-''^^'iUt=iPi)]- Using the 
identity A~'^ exp{B) A = exp{A~^B A) we get 



Uitc) = n ^''""^ 



(2.3) 
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The operators Hk are the Hamihonians in the so-caUed 
"togghng frame" jll| 

k 

Uk ■■= Qk^HQk, with Qk =l[Pi (2.4) 

1=1 



By expanding the exponential in (2.3) and collecting 
terms of equal order in TikTk the time evolution oper- 
ator can be written as 



[/(i,) = e-(«+«'^'+«<^'+-K, 



(2.5) 



where the operators 7i, 'H^^\ . . . are the average Hamil- 
tonians of increasing order. In the remainder of this pa- 
per we will only use pulse sequences consisting of equally 
spaced pulses (r^ = t) . For such a sequence with n pulses 
( "nr-sequence" ) the leading order terms of the AH ex- 
pansion are given by 
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(For a more detailed derivation refer to App. B of p2[.) 
The time evolution can therefore be described approxi- 
mately by the average of the toggling frame Hamiltoni- 
ans, 7Y, in the following simply referred to as "the average 
Hamiltonian" . This holds exactly for vanishing commu- 
tators between all Tik- 

For a given time period fc, the pulse sequence will lead 
to the transformation Qk — (S'/i^i Qk 
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.^.=1-.^ , where Q^-' is 
the transformation on qubit /i. The applied pulses, and 
hence the operators Qj^ , are unitary transformations. 
Any unitary transformation in SU(2) can be represented 
by ordinary 3-dimensional rotation matrices R = (Rij) 
via the relation 






(2.8) 



so that any pulse sequence can be characterized by a set 
of matrices R}/* . 

The resulting average Hamiltonian can be writt en i n 
the same form as the original system Hamiltonian (2.1); 
it is straightforward to verify from (E^) and (|2.8D that 



the vector representing the Zeeman part for qubit fi is 
given by 
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and the 3x3 tensor representing the coupling part for a 
given pair /i, i/ by 



H^'' = 
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(2.10) 



Our goal is to effectively remove all interactions by ap- 
plying a proper pulse sequence. The more general task 
of generating an evolution corresponding to selected cou- 
plings being turned on or off requires only minor mod- 
ifications. Specifically, we want to achieve "first order 
decoupling" , meaning that the average Hamiltonian van- 
ishes. For now higher order terms in the AH expansion 
will be neglected. Also, the Zeeman interactions shall be 
excluded for the time being. 



III. DECOUPLING FOR 2-QUBIT NETWORK 

We consider at first the case of two qubits. Ignoring 
Zeeman interactions we are left with a single interaction 
term represented by the tensor H^^ =: H. First order 
decoupling is then achieved by a pulse sequence resulting 
in H = 0. 

Two principally different kinds of pulses can be used, 
spin-selective and nonselective ones. From ( [2.10 ) the fol- 
lowing theorem on the selectivity of the required decou- 
pling pulses can be derived: 



Theorem 1 A coupling represented by a tensor H with 
TtH ^ cannot be decoupled with nonselective pulses. 



Proof: For nonselective pulses (Rj. 



(2) 

\j. ~ Rfe ) the toggling 
frame coupling tensor is an orthogonal transformation of 
the original coupling tensor, therefore Tr H^ = Tr H = 
TrH. For decoupling one needs H = 0, and hence Tr H 
must vanis h. (T his also holds for sequences with different 
Tk, where (2.1L) is then given by a weighted average.) D 



In the case of, e.g., direct dipolar couplings (with 
H oc diag(l,l,— 2) in the high-field approximation) this 
condition is fulfilled and nonselective sequences like the 
well-known WHH-4 [Q can be applied. On the other 
hand, the widely discussed strong or weak scalar cou- 
pling requires selective pulses. 

We now give a specific first order decoupling sequence 
for arbitrary couplings using spin-selective tt pulses with 
respect to the coordinate axes x,y,z. These pulses are 
given by the operators X = ^-^'""'^1'^ = \a^^ and analo- 
gously for Y and Z . The corresponding rotation matrices 



X = R, -diag(+l,-l,-l), 
r = Rj, = diag(-l,+l,-l), 
Z=R, = diag(-l,-l,+l), 



(3.1) 



with the multiphcation table shown in Table |. The de- 
coupling is based on the following property 



aol + aiX + a2Y + asZ = ^ 



(3.2) 



where / denotes the identity matrix. A weighted sum of 
the matrices /, X, Y, Z vanishes if and only if the weight 
factors are equal. 

Following B a pulse sequence will be written as a ma- 
trix of transformations /, X, Y, and Z specifying the 
corresponding Q or R on a specific qubit for a certain 
time interval of duration r. The n time intervals corre- 
spond to the columns (from left to right) and the differ- 
ent qubits fj, — 1,2, ... ,N to different rows. We may thus 
speak of an n-dimensional row vector R-^' , the elements 
of which are the ordered sequence of transformation ma- 
trices for given /j,, RJ/* , k = 1,2, ... ,N. This represen- 
tation gives more direct insight into what is happening 
to the coupling Hamiltonian than the notation involving 
the actually applied pulses. 

Now consider the sequence 






I I I I 
I X Z Y 



(3.3) 



which can be implemented by the following train of 
TT pulses on qubit 2 [cf. (|2.2|)] 



(3.4) 



(T,x(2),r,y(2),T,x(2),r,-y(2)). 



(The final y-pulse is applied in order to make the se- 
quence cyclic.) The resulting average Hamiltonian is 
given by 

H^-HYr'^^^ =^-H{I + X + Y + Z), (3.5) 



which vanishes due to (3.2) for any coupling tensor H. 
This sequence thus effectively removes an arbitrary cou- 
pling between two qubits. The Hamiltonian does not 
even have to be known, because the same sequence works 
for any interaction. Hence it follows: 

Theorem 2 Any interaction between two spins can be 
decoupled to first order using a sequence of spin-selective 
pulses. 



For diagonal interaction models (2.10) becomes 



-HERi^^Ri'^=:-HVR 



(12) 
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because all matrices are diagonal (for the n pulses under 
investigation). Then it does not matter how the trans- 
formations are "distributed" between the two qubits, be- 
cause only the product needs to be considered. Thus 



alternative sequences are possible, which will be particu- 
larly interesting for the N > 2 case discussed later. The 
sequence 
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X 
Y 



(3.7) 



lea ds to the same transformations of the coupling tensor 
as (3^) and thus also decouples the interaction. Fur- 
thermore, it has the additional feature of removing the 
Zeeman terms, for the typical case that they only con- 
sist of a z-component, tVzd'z. The corresponding pulses 
applied both on qubit 1 and 2 are given by 



{t,X 
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The sequential order of the transformations clearly has 
no influence on H and thus on the first order decoupling. 
For the above mentioned sequences it is chosen so that 
no z-pulses are needed, which usually cannot be imple- 
mented in conventional NMR. The higher order terms, 
however, are affected by the order. 



IV. DECOUPLING FOR Ai'-QUBIT NETWORK 

We now turn to a system consisting of N qubits. It 
often suffices to take only nearest and, possibly, next 
nearest neighbors into account. We will come back to 
this issue later. For now we discuss the most general 
case of a completely coupled system, i. e. each qubit is 
coupled to all the other ones. From (2.10) it is evident 
that each pair interaction H^'' can be treated separately, 
essentially leading back to the N ^ 2 case. So, we can 
use the sequence for A'^ = 2 as a starting point. An ut- 
decoupling sequence for A^ qubits will be represented by 
an A X n matrix. Decoupling means elimination of all 
H'^'^. The task is to find a matrix that yields a valid 
N — 2 decoupling sequence for any pairing of rows. 



A. Simple Scheme 

The simplest approach is to recursively nest the se- 
quence ( ^.31 ) in analogy to the method in H]. When a 
qubit is added to a valid sequence, it will be exposed to 
a complete unit cycle during each time interval of the 
initial sequence. For 3 qubits we thus get: 



I I I I I I I I I I I I I I I I 

II I IXXXXZZZZYYYY 
I X ZY Y Z X I I X ZY Y Z X I 



(4.1) 



While the additional qubit will undergo one 4r unit cy- 
cle, the transformations on the other qubits are constant, 
and it will therefore decouple. However, this scheme is 
not efficient in terms of the time intervals r and pulses 
needed, which both scale as 4^~^ because the sequence 



has to be prolonged by a factor 4 for each added qubit. 
The advantage is that at each time interval only one pulse 
has to be applied. This is why every second unit cycle is 
reversed. Furthermore, the almost realized mirror sym- 
metry {Tik — Hn+i-k) is expected to reduce the higher 
order correction terms: for a symmetric sequence, all cor- 
rection terms of odd order vanish f|12| . 



B. Efficient Scheme Using Orthogonal Arrays 

Much more efficient schemes can be constructed using 
the notion of orthogonal arrays (OA), which have been 
extensively studied in the context of combinatorics, de- 
sign theory, and error-correcting codes. A comprehensive 
treatment can be found in p4| . First we want to give a 
definition of orthogonal arrays and then demonstrate how 
they are related to the decoupling problem. 

Let S* be a set of s symbols. Then an orthogonal array 
is defined by 

Definition 1 A k x As* array A with entries from S is 
said to be an orthogonal array with s levels, strength t, 
and index X {0 < t < k) if every t x As* subarray of 
A contains each t-tuple based on S exactly A times as a 
column. It will be denoted by OA(Xs* ,k, s,t). 

For our purposes we will only regard the case t — 2 and 
s — 4. Usually the symbols are chosen as the numbers 0, 
1, 2, 3. We will interpret them as the n rotation matrices 
including the identity matrix, so that S = {I, X,Y, Z}. 
As an example regard the following QA(16, 5, 4, 2) 



/// IIXXXXYYYYZZZZ\ 
IXYZIXYZIXYZIXYZ^ 
I X Y Z X I Z Y Y Z I X Z Y X I 
IXYZYZIXZYXIXI Z Y 

\lXYZZYXIXI Z Y Y Z I X J 



(4.2) 



If you pick any two rows (e. g. row 2 and 5), each of the 
s* — 16 possible 2-dimensional column vectors appears 
exactly once, hence A = 1. Definition nl implies fur- 
ther that every row contains each symbol exactly As = 4 
times. As it turns out, these properties are sufficient for 
the array to be a valid decoupling sequence: 

Theorem 3 Any QA(16A, fc,4, 2) can be used to decouple 
fc + 1 qubits within n ~ 16A time intervals. 

Proof: Pick two rows /i and v of the OA and consi der th e 
corresponding average Hamiltonian according to (2.1C), 



)M 



- 16A. Each of the s* = 

(/,/), (/,X),...,(Z,Z) 



16 combinations (R^'' ^R^'^ ) 
occurs exactly A times. We thus get H^'' — j^{I+X+Y+ 
Z) H'"^ {I + X + Y + Z), which vanishes for any H^". We 
can further add an additional (k + l)th row to the array, 
consisting only of I 's. As every row of the OA contains 
an equal number Xs of each symbol, the AH involving the 
{k + l)th qubit is H''''=+i = iH^-'^+i {I + X + Y + Z)=0. 
Thus, the OA can be used to decouple fc + 1 qubits. D 



Using large enough OAs (fc -I- 1 > N), decoupling se- 
quences for an arbitrary number N of qubits can thus be 
constructed. We note that simultaneous pulses on dif- 
ferent qubits are necessary for this scheme to work. The 
essential question is for which sets of parameters A, fc, s, t 
an OA exists at all. Specifically, we want to find arrays 
with fc as large as possible, in order to decouple many 
qubits with the least effort. 

Bose and Bush |13] have given a recursive construction 
scheme for arrays of strength 2, providing OAs with large 
values fc. From their Theorem 4 follows the existence of 
all Q4(n, fc,4,2) with 



i(n-l), u = 0,2,4,. 
i(n — 5), u — 1, 3, 5, . 



(4.3) 



An example of the underlying construction scheme will be 
given in the appendix. This guarantees the existence of 
the arrays 04(16, 5, 4, 2), 04(32, 9, 4, 2), OA(64, 21, 4, 2), 
04(128, 41, 4, 2), and so on. [Note that this construction 
does not yield all possible O As, t here exists, for example, 
an QA(48, 13,4, 2).] From (4.3) we can derive a lower 
bound on the efficiency of the scheme: An [n — 2"+^)t- 
sequence can be used to decouple up to A^ = i(n ± 2) 
qubits (the sign depending on u being even or odd). It 
follows that N qubits {N > 3) can be decoupled by a 
sequence consisting of ri = cN — 2 time intervals, where 
3 < c < 6. The number of pulses that are needed for 
the sequence is less than N{cN — 2). Hence it is possible 
to decouple (to first order) any interaction with linearly 
increasing effort in time: the scheme is efficient. 



C. Efficient Scheme for Diagonal Couplings 

For diagonal couplings there exists an even more effi- 
cient scheme. When the transformations on two qubits 
/i a nd V are given by the row vectors W^' and W^' 
[cf. (|3.3|)], the average Hamiltonian is, in analogy to (p^). 



V\^"' = -H^^^ R^'^^R^^'^^ =: -H^'' (il^^^- fi^")) . (4.4) 
"" fc=i ^ 

It vanishes whenever the "scalar" product R-^' ■ i?*-"' 
(which is really tensor valued) equals zero. This leads 
with (^^) to the following decoupling criterion: 

Theorem 4 A pulse sequence represented by a matrix of 
N rows can be used to decouple N qubits with diagonal 
coupling between any qubit pair if the scalar products be- 
tween different row vectors vanish. This is equivalent to 
the condition that the vector of the element-wise prod- 
ucts of two rows contains each of the elements I, X, Y, Z 
equally often. 



The sequences already discussed, like e.g. (4^), ful- 
fill this condition, but there exist suitable matrices with 
a smaller number of columns, leading to even shorter 



sequences. These matrices are closely related to the so- 
called difference schemes or matrices [Q, which play an 
important role in the construction of the OAs discussed 
in the previous section. Difference schemes are defined 
over a finite Abelian group of s elements with a binary 
operation +, (.4,+). 

Definition 2 An r x c array D{r, c, s) with entries from 
A is called a difference scheme based on {A, +), if for all 
fx and V with 1 < ^i^v < r, ^ ^ v, the vector difference 
between the fith and vth rows contains every element of 
A equally often. 

Relevant here are difference schemes over the additive 
group of the Galois field GF{s) (cf. App. A of 0) with 
s = 4 elements: 

Theorem 5 A D{r, c, 4) over (GF(4), +) is equivalent to 
a CT-decoupling sequence for N ~ r qubits. 

Proof: The group of tt rotations, {{I, X.Y, Z},^), is iso- 
morphic to (GF'(4),+), which is evident from TableU. So 
the elements I , X , Y , Z can be identified with 0, 1, 2, 3, 
and the matrix product with the operation +. Each ele- 
ment is its own inverse, so that the difference of two rows 
equals the sum. Then the defining property for difference 
schemes is equivalent to the condition of Theorem 0. D 

For difference schemes one has r <c ll^l, so that 
at most N — n qubits can be decoupled with an nr- 
sequence. An r x r array D[r,r,s) is also called gen- 
eralized Hadamard matrix. 

Multiplying a column with a specific element leaves 
the scalar products invariant, so that the matrices can be 
brought into a normalized form, with only /'s in the first 
row. The other rows then have to consist of permutations 
of equal numbers of /, X, Y, Z, so that from the outset 
only sequences with n = 4u (m £ N) time intervals are 
applicable. By trying all combinations we have found the 
following matrices for decoupling of up to 4 or 8 qubits, 
respectively. 




M« = 



Ma 



( I I I I I I I I \ 

I I X X Y Y Z Z 

I X Y Z I X Y Z 

I X Z Y Y Z X I 

I Y I Y Z X Z X 

I Y X Z X Z I Y 

I Z Y X Z I X Y 

\l Z Z I X Y Y X ) 



(4.5) 



(4.6) 



the also existing M12. One is led to suspect that there 
may exist matrices M„ for all n = 4u, m g N, which, 
however, is still an open question. But at least for all 
n = 2", M > 2, difference schemes can be constructed 
using the properties of Galois fields [ISl. (For a more 
detailed discussion on the existence of difference schemes 
see chap. 6 of [|lj|.) A simpler way to find matrices M„ 
can be based on the following theorem: 

Theorem 6 Given two decoupling matrices A/„j, Mn^, 
their direct product M„-^n2 = -^ni "X" M„2 is again a de- 
coupling matrix. 

Proof: The (712(1 — 1) -t- k)th row {i — l,...,ni;k — 
1, . . . ,712) of the matrix Mn^n^, ihk), derives from the 



ith row vector of Aim, 
k, as (i, k) — (iiki . 



and the kth row vector of Mn 



iik„ 



Znjfci 



7ni "-712 J 



The scalar product between two rows (i, k) and (j, I) is 

{hk)-{j,l) ^iiji{k-l)-\ VinJnA^-l) = {i-3){k-l)- 

Because all rows in Af„j and Mn^ are orthogonal to each 
other, (i, k)-{j, I) can only be different from zero for j = i, 
I = k. Therefore all rows o/M„j„2 '^'"^ orthogonal. D 

Now using the matrices Af4 and Mg given above one 
can construct decoupling matrices M2U for any u > 2. 
Consequently, for N qubits there exists a {cN)t- 
decoupling sequence with at most cN'^ pulses, where 
1 < c < 2. So, for diagonal couplings the efficiency is 
further improved compared to the general case of the 
last section. 



D. Zeeman terms and partially coupled systems 

Removing all Zeeman terms requires only a small mod- 
ification. From (2.9) it follows that the Zeeman terms of 



We will denote such n x n decoupling matrices corre- 
sponding to a D{n,n,4) by M„. We do not show here 



a specific qubit 11 are removed if the transformations on 

that qubit sum up to zero, J2k=i ^k ~ ^- "^^^^ corre- 
sponds to a vanishing row sum in the matrix of trans- 
formations. For the given sequences this is fulfilled for 
all rows except the first. Thus only the Zeeman interac- 
tion of qubit 1 is not removed. So, all Zeeman terms can 
be suppressed by using a matrix for iV -I- 1 qubits and 
omitting the first row. 

Similarly, leaving selected couplings intact rather than 
removing all interaction terms takes only small alter- 
ations. To achieve, for example, a time evolution cor- 
responding to one single coupling between qubits 1 and 
/I, it is merely necessary to use the same transformations 
on qubit /i as are used on qubit 1. 

In reality the couplings are usually well localized and 
it is sufficient to take only close neighbors into account. 
Then much simpler sequences can be used. The situation 
essentially can be regarded as a coloring problem M . For 
nearest-neighbor couplings the problem reduces to the 
N = 2 case, so that (3.3) or (3.7) can be used, subjecting 



all qubits with odd number to the same transformations 
and likewise for the even qubits. 



V. HIGHER ORDER TERMS 



with 



Up to now we have only considered first order decou- 
phng, neglecting the higher order terms in the AH expan- 
sion. These terms can be suppressed further by repeating 
a given sequence with correspondingly shorter time inter- 
vals. 

Let us regard the sequence we get from repeating m 
times a given nr base sequence during a given cycle time 
tc- It then consists of mn time intervals of duration 
tc/{mn) each. The average Hamiltonian is the same as 
for the base sequence and v anis hes, while for the first 
correction term we get from (^.7]) 



n 



(1) _ 



itc 



mn 



j>k 



[Hj , Hk 



(5.1) 



By using the new indices s = l,...,m for the 
different subcycles and k = 1 , . . . , n for the time 
interval within a given subcycle the sum can be 
written as Er>"fe [^j.^fc] = Er=s' Ei>fe [^j . ^fe] + 
J27>s' YTj.k [^h'^k] =mJ2'i>k [^j>^fc]: implying 



H(i) = - 



2mn' 



Y.^n,,nk]. 



(5.2) 



i>k 



The first correction term is thus reduced by the factor 
1/m if the base sequence is repeated m times. In prin- 
ciple, by using arbitrarily fast pulses, perfect decoupling 
could be realized. 
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(a'-hi) 
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a['^\ 



G 



(A'+l)^;i(A'+2) 



<yi 



't] 1 '■^"] ' -'tjk-"i '^"] ^"k ' 

As expected, the correction term vanishes for weak cou- 
pling, where only a^ terms are involved (Ising model, 
Jx = Jy = 0) and the ti me e volution is given exactly by 
the average Hamiltonian ( |2.6| ). We now define the fidelity 
of a given initial state \'ip) after the time tc 



F{U)^ {'^\U[tc)\'^) 



(6.3) 



As an example we restrict ourselves to strong scalar cou- 
pling (Heisenberg model, Jx = Jy = Jz — J) and equal 
chemical shifts (wz = w); taking as the initial state the 
state with qubit 1 excited and all other qubits in the 
ground state, |i/)) — |100 . . .0), the fidelity for small tc is 
approximately 



ntc) 



1- 
1- 






[{J^ + oo^){N - 1) 



iV = 2, 

2 J2] , iV > 3. 
(6.4) 



In Fig. |i| the undisturbed evolution is plotted for A^ = 4 
qubits together with the curves for two decoupling se- 
quences with different pulse frequency m. It can be seen 
how increasing m improves the fidelity, i. e. the suppres- 
sion of the higher order terms. 



VI. EXAMPLE: SPIN CHAIN 

We will now study a simple example, a spin chain 
consisting of N spins with a diagonal coupling between 
nearest neighbors. The corresponding Hamiltonian is 
[cf. (O)] 
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N 
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c.l^)al^) 



^(J.4^)®a(' 



s.(m+i) 



J,a(^)( 



^(m+1) 



JzO- 



^(A'). 
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(6.1) 



For this interaction model, the N — 2 sequence for di- 
agonal couphngs (3.7,3.8) can be used (thus n — 4), ap- 
plying the a;-pulses on all qubits with odd number and 
the 2/-pulses on all qubits with even number. The lead- 
ing correction term Ti,^^' does not vanis h, in general: for 
the respective sum of commutators in (2.7) one gets by 
induction 




FIG. 1. Decoupling for A^ = 4 qubits, lj = 0. Shown is the 

fidelity F{tc) for free evolution ( ), decoupling with m = 1 

( ) and the corresponding approximation F = 1 — ( Jt)'*/4 

( ), decoupfing with m = 4 ( ) and the approximation 

F=l~{Jtf/64: ( ). 
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+ Y.{^f^-^^^i'i'+^^9bA't>) \, (6-2) 



VII. DECOUPLING WITH FAULTY GATES 

Up to now we have only considered sequences consist- 
ing of perfect pulses, implementing exactly a tt rotation. 
This assumption is, of course, not tenable in any real 



physical system. Even for ideal pulses the decoupling is 
generally not perfect due to higher order terms, but as 
stated above, these terms can be suppressed by increasing 
the pulse frequency m. Yet with faulty gates this strat- 
egy is not necessarily successful, because more errors are 
then introduced into the system from the outside. 

In order to see what happens we study a very simple 
model, two spins that interact via scalar coupling. We 
again use sequence ( |3.8| ) but now we assume faulty x- 
pulses, where the rotation angle differs from tt by a small 
deviation 6, so that 



X^^\S) 



i ai^' cos ■ 



1 sin ■ 



(7.1) 



where S is randomly distributed with the probability dis- 
tribution wis, a) = (V27^CT)"^e"''''/2'^^ The standard 
deviation a measures the scatter of the pulses. 
As a first case we consider weak coupling, with 



W = -c.«a« -4'^a(2) + Ja«al2). 



(7.2) 



Of course, in this case one would not use the above strat- 
egy because the higher order terms vanish anyway, and 
increasing the number of applied pulses will, with cer- 
tainty, lead to a degraded performance. By induction 
one can show analytically that after application of m de- 
coupling cycles the fidelity for one of the product basis 
states (|00), |01), |10), |11)) as initial state is given by 



Faim) = - (1 + e 



(7.3) 



For increasing a, i. e. faultier gates, the fidelity decreases 
as expected. 

For the case of strong coupling. 



n 



(aW+af)+J^(i).^(2)^ (7 4) 



we have resorted to a numerical calculation of the fidelity 
after the time tc for the initial state |01), for different a 
and numbers of repetitions m of base cycles. For each set 
of parameters we calculated the mean over 1000 realiza- 
tions with randomly picked deviations 5. Our results are 
shown in Fig. g. It turns out that the resulting behav- 
ior can be approximated pretty well by the product of 
the function for a = 0, approximated by (|6.4|) , with the 
"damping" resulting for weak interaction (|7.3|), yielding 
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1- 
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400 m2 



(7.5) 



As one can see, we have indeed two competing effects: for 
small m and a the suppression of the higher order correc- 
tions is dominant, whereas for larger m and a the errors 
introduced by the faulty gates dominate. In particular, 
there exists for a given a an ideal number of iterations m 
which yields the maximum fidelity. The corresponding 
^max decreases thereby for increasing error. This means 



that for faulty gates it is not possible anymore to achieve 
perfect decoupling. This conclusion gets increasingly sig- 
nificant for larger N: then i^max is further reduced by 
the effect of the additional faulty gates and the rising 
influence of higher order terms. 
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FIG. 2. Decoupling for N — 2 qubits with faulty gates, 
strong J-coupling. Shown is the fidelity F^ at tc = 1/J, 
J = lOoj. The curves correspond to a = (•), 0.01 (x), 
0.02 (A), 0.03 (•), 0.04 (a), 0.05 (♦). The approxima- 



tion function is F^{m) = \{l + e-"'' )(1 - j^ 
with correlation coefficient 0.999996 for a = 0. 



+ 



1369 m* ''' 



VIII. SUMMARY AND CONCLUSIONS 

We have shown that it is possible to generalize refo- 
cusing methods to systems of spins coupled by arbitrary 
two-spin interactions. Thereto spin-selective operations 
are necessary. We have given a concrete sequence for first 
order decoupling for N = 2, and schemes for constructing 
sequences for A^ > 2. These schemes are efficient regard- 
ing the needed time, which scales linearly with the size 
of the system, and the amount of pulses, which scales at 
most quadratically. 

However, in the general case there are higher order 
terms that must be taken into account. These can be sup- 
pressed by increasing the pulse frequency, assuming ideal 
gate operations. Perfect decoupling for long times could 
only be achieved by infinitely fast decoupling pulses. 

These results hold for the most demanding case of 
an arbitrary coupling model. Yet for special interaction 
types, much more efficient schemes are possible. Direct 
dipolar couplings, for example, can be decoupled using 
nonselective pulse sequences, that work independent of 
the system size. For weak scalar coupling all higher or- 
der correction terms vanish, so that the decoupling is 
perfect. The sequences can also be considerably simpli- 
fied for only partly coupled systems, like a spin chain 
with nearest neighbor couplings only. 

In a real physical system pulses never are ideal. For 
faulty gate operations perfect decoupling is no longer pos- 
sible, because of the inaccuracies introduced into the sys- 
tem by the procedure itself. 



The present results are not restricted to the special 
task of removing all interactions in a spin system. Minor 
modifications allow specific coupling terms to be turned 
on or off. Combining on- and off-periods one can even 
rescale the various parameters within the original Hamil- 
tonian. Generally, a wide range of effective Hamiltonians 
can thus be invoked. In this sense, complete decoupling 
is a special case of a quantum simulation [|f 7| , where a 
system is modeled that "does nothing" , i. e. any initial 
state remains unchanged. Such a system is, in general, 
quite challenging to realize and may thus be considered a 
serious testing scenario for quantum dynamical control. 



ACKNOWLEDGEMENTS 

We thank J. Gemmer, I. Kim, A. Otte, F. Tonner, and 
T. Wahl for fruitful discussions. 



APPENDIX: RECURSIVE CONSTRUCTION OF 
ORTHOGONAL ARRAYS 

This construction is due to Bose and Bush flia. Ex- 



amples can also be found in 14 



1. Construct an Q4(16A,4A,4, 2), A = 2", u e N, as 
follows: Take an M^x [a i:'(4A,4A,4)]; then the OA 
Aq is given by 

Ao = M4A «) ( / X Y Z). 

2. If possible, construct an OA A\ with Ai — A/4, ac- 
cording to 1. Repeat each column 4 times, yielding 



K = M 



I I I I 



If possible, repeat this step with an OA A2 where 
A2 = A/4^, and form an A2 by repeating each col- 
umn 4^ times, and so on. Append the matrices A[, 
A'2, ... to Aq. 

3. Add a row consisting of 4A /'s, 4A X's, and so on. 

We now have an Q4(16A, k, 4, 2) with A = 2", u G N, and 
fc = 4A + 4Ai + • • • + 1 rows, or 



A(4= 



1) 



ic-l 



m/2, u = 0,2,. 

{u-l)/2, u = l,3,. 



The OA (4.2) has been constructed in this way, using 
M4. A Ai = 1/4 does not exist, so step 2 was omitted, 
resulting in fc = 4 + 1 = 5. An QA(32, 9,4, 2) can be 
constructed using Mg. For an QA(64, 21,4, 2) we con- 
struct ^0 using M16 = M4 (g) M4 and Ai using M4, so 
fc = 16 -F 4 -h 1 = 21. 
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TABLE I. Multiplication table for the matrices /, X, Y, Z 
and addition table for GF(4) 
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